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Abstract. 

In two recent INTEGRAL papers, Lutovinov et al. (2003) and Blay et al. (2004) report a timing and spectral 
analysis of the transient Be/X-ray pulsar SAX J2103.5+4545 at high energies (5-200 keV). In this work we present 
for the first time a study of the pulse profile at energies above 20 keV using INTEGRAL data. The spin-pulse 
profile shows a prominent (with a duty cycle of 14%) and broad (with a FWHM of ~ 51 s) peak and a secondary 
peak which becomes more evident above 20 keV. The pulsed fraction increases with energy from ~ 45% at 5-40 
keV to ~ 80% at 40-80 keV. The morphology of the pulse profile also changes as a function of energy, consistent 
with variations in the spectral components that are visible in the pulse phase resolved spectra. A study of the 
double peaked profile shows that the difference in the two peaks can be modeled by a different scattering fraction 
between the radiation from the two magnetic poles. 
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1. Introduction 

After the discovery by Hulleman et al. (1998) and sug¬ 
gestion that it is a Be X-ray binary, SAX J2103.5+4545 
was conclusive confirmed as such by Reig et al. (2004) 
who extimate a distance of ~ 6.5 kpc. The system con¬ 
sists of an accreting neutron star (NS) and a BOV main- 
sequence donor star which shows Balrner emission lines 
(Reig et al. 2004). The relatively short orbital period, 12.7 
days (Baykal et al. 2000), and the spin period of ~ 358 s 
(Hulleman et al. 1998) make SAX J2103.5+4545 a pecu¬ 
liar Be/X-ray transient system. The source does not follow 
the Pp U i S e versus P or b relation followed by all other Be 
X-ray binaries as initially established by Corbet (1986). 
The source has been extensively studied in outburst by 
Hulleman et al. (1998) with BeppoSAX , Baykal et al. 
(2000, 2002) with RXTE, inam et al. (2004) with RXTE 
and XMM-Newton, Reig et al. (2004) at optical and in¬ 
frared wavelengths, and by Lutovinov et al. (2003), Sidoli 
et al. (2004) and Blay et al. (2004) with INTEGRAL. 

We present here for the first time spin phase resolved 
measurements at energies above 20 keV, based on hard X- 
ray INTEGRAL observations of SAX J2103.5+4545 dur¬ 
ing a bright outburst. 
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2. Results 


The data were obtained during the AOl INTEGRAL ob¬ 
servation of the Cygnus region, performed from May 2 to 
May 3 2003 (52761.32-52762.43 MJD). We use data from 
the coded mask imager IBIS/ISGRI (Ubertini et al. 2003; 
Lebrun et al. 2003) with an exposure of 75.06 ks, and from 
the JEM-X monitor (Lund et al. 2003) with an exposure 
time of 33.34 ks. The data reduction was performed us¬ 
ing the standard software and response files provided in 
OSA 3.0 ( Goldwurm et al. 2003 1. We extracted 32 chan¬ 
nel spectra for IBIS/ISGRI and 256 channel spectra for 
JEM-X. The spectra were analyzed with the XSPEC V. 
11.3 software package (lArnaud 19961) . The background 
subtracted ISGRI light curve for the timing analysis was 
reduced using dedicated software ( iiJight_extract , ver¬ 
sion to be included in OSA 5.0). The ISGRI light curve 
per energy band was obtained through a simultaneous fit, 
for every 20 s long time bin, of all sources detected in 
the overall image and the background using the standard 
least squares method in the detector domain. For JEM-X 
the background subtracted light curve was obtained using 
standard software (see JEM-X Analysis User Manual). 
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2.1. X-ray Timing 

We searched for coherent pulsations of the source in the 
20-40 keV energy band where we have the best statis¬ 
tics. Using the 20 s binned light curve of the source in the 
20-40 keV energy band, we computed a Power Density 
Spectrum (PDS) in the frequency range between 5xl0 -5 
and 0.025 Hz from Fast Fourier Transforms. In the result¬ 
ing PDS an evident signal is present at v = 2.8296 x 10 -3 
Hz, and first and second harmonics are also visible. The 
highest peak in the power spectrum, the first harmonic, 
corresponds to a nominal period of 353.41 s. The accurate 
period was found using the phase-delay fitting method (see 
e.g. JNagase 1989[ > resulting in P = 353.45 ± 0.16 s, with 
an upper limit to the pulse period derivative of |P| < 2.2 s 
yr . In this section the errors are given at ler confidence 
level for a single parameter. The measured pulse period 
is in agreement with the expected value using the pulse 
period of 354.794 s (52761.87 MJD) and the spin-up rate 
of ~ 7.4 x 10 -13 Hz s _1 reported by inarn et al. (2004). 



0 0.5 1 1.5 2 


Pulse phase 

Fig. 1. JEM-X and IBIS/ISGRI light curves of SAX 
J2103.5+4545 folded at the NS spin period. From top 
to bottom, in the energy range 5-20 keV (JEM-X), 20- 
40 keV and 40-80 keV (IBIS/ISGRI), respectively. Each 
pulse profile is repeated once for clarity. 

Fig. □ shows the folded light curve at the best pe¬ 
riod. These profiles show a very prominent peak (~ 51 s 
FWHM) around phase ~ 0.8. The pulse fraction (PF), de¬ 
fined here as PF = (Imax Imin)/(.Umax + Imin), increases 
with energy, going from 47.5+0.3% in the 5-20 keV energy 
band to 45.4+6.3% at 20-40 keV and to 80.2+8.9% at 
40-80 keV. In the energy band 80-150 keV the low signal- 
to-noise ratio does not allow a pulse fraction to be deter¬ 
mined. The X-ray pulse profile of SAX J2103.5+4545 also 
shows an energy dependent secondary peak around phase 
ss 0.2, which becomes more evident at energies above 20 
keV. Finally, note that the hick-ups in the top panel at 
phase 0.6 seem to be consistent with those observed in 
the pulse profile of this source by RXTE/XMM-Newton 


(e.g., Inam et al. 2004). Note that these hick-ups are also 
present in the 40-80 keV pulse profile (bottom panel), 
while are not significantly present in the 20-40 keV pulse 
profile (middle panel). 

2.2. Spectral analysis and pulse phase resolved spectra 



Energy (keV) 


Fig. 2. Ratio of the “main” to “secondary” peak of the 
pulse profile as a function of energy from 5 to 200 keV. 
The dotted line is from the best fit two polar cap diffusion 
model (see §3). 

We extracted the averaged over the whole observation 
spectrum of the source. In order to investigate the pulse 
phase dependence of the high energy emission from this 
source, we selected two phase intervals in the pulse profile 
of Fig. [I] one corresponding to the main peak (0.5-1.0) 
and the other corresponding to the secondary peak (0- 
0.5). From each of these phase intervals we extracted JEM- 
X and IBIS/ISGRI spectra. The spectral analysis has been 
carried out in the 5-20 keV energy range for JEM-X and in 
the 20-200 keV energy range for IBIS/ISGRI; the JEM-X 
spectra were rebinned to 28 channels to limit the oversam¬ 
pling of the instrumen resolution to a factor of 3. 

The two phase-resolved spectra look different as is ev¬ 
ident from their ratio plotted as a function of energy in 
Fig. El The ratio gradually increases from 5 to 30 keV and 
then decreases at higher energies. Note that the Direct 
Continuum (DC) level is not subtracted from these spec¬ 
tra; this might introduce a variation of the spectral ratio 
simply due to the fact that the DC level changes with 
energy. In principle we could estimate the DC level from 
the spectrum at the minimum of the folded light curve, 
although we do not know if all the emission at the min¬ 
imum of the light curve is due to a DC level, or instead 
is due to unpulsed emission from the polar caps (this is 
possible for instance if parts of the polar caps are always 
visible). However, the statistics of these data is not good 
enough to extract a meaningful spectrum at the minimum 
of the folded light curve. Note that we took into account 
the presence of a DC level in the diffusion model that we 
used to fit this spectral ratio (see §3). 
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Table 1. Best fit parameters of the absorbed Comptt model for the spectra of SAX J2103. 5+4545 together with 
90% confidence errors for a single parameter. 


Pulse phase 

0.0 - 0.5 

0.5 - 1.0 

All 

All 

A h (10^cm _ “, fixed) 

3.45 

3.45 

3.45 

3.45 

kT 0 (keV) 

1.5 (f) 

1.5 (f) 

1.50 ± 0.17 

1.5 (f) 

kT e (keV) 

13.0 (f) 

13.0 (f) 

13.68 ± 1.85 

13.0 (f) 

T 

2.08 ± 0.17 

2.5 ± 0.13 

2.37 ± 0.43 

2.53 ± 0.14 

Xred 

1.2 (49 d.o.f.) 

1.2 (49 d.o.f.) 

1.09 (47 d.o.f.) 

1.06 (49 d.o.f.) 

L%-20QkeV (10 36 ergs 3 ) 

2.95 

4.64 

4.17 

4.21 


* Assuming a distance of 6.5 kpc. 
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Fig. 3. The 5-200 keV unfolded spectrum of SAX 
J2103.5+4545 along with the best fit COMPTT model. 


To fit the total spectrum, we first used an absorbed 
cut-off power law model (xi e< i = 1-28, 48 d.o.f), were the 
photon index was found to be T ~ 1.45 with an cut-off at 
E cu t ~ 30 keV. A better fit was found using an absorbed 
power law multiplied by a high-energy exponential cutoff 
model (Xred = 1-37, 47 d.o.f.), where the parameters found 
were similar to the ones reported in Blay et al. (2004). 
We then substituted this phenomenological model with 
a physically motivated thermal Comptonization model 
COMPTT, using the slab geometry (Titarchuk 1994), which 
gives a yfed °f 1-09 (47 d.o.f.). Usually the power law 
plus high-energy cutoff model is a good approximation 
of the spectra of high-magnetic field X-ray pulsars, due 
to the fact that the strong magnetic field modifies the 
Compton cross section giving a sharp cutoff in the spec¬ 
trum. However, in our case, the COMPTT model provides 
a slightly better description of the data with respect to 
the “empirical” power law multiplied by a high-energy ex¬ 
ponential cutoff on the entire dataset at 81% confidence 
level (estimated by means of an F-test). We therefore pre¬ 
fer to report here the best-fit parameters of the COMPTT 
model. We interpret the parameters of the COMPTT model 
in the following way: at the bottom of the accretion col¬ 
umn, i.e. at the hot spot on the NS surface, a blackbody 
emission, with a temperature fcTo, contributes the soft 
component of the spectrum, and provi des also t he seed 
photons for the Comptonization lllnam et ah 2004). These 
are Comptonized in the hot plasma present in the accre¬ 
tion column at a temperature kT e and optical depth r. 


In fitting the pulse phase resolved spectra, we found 
that the spectral parameters of the COMPTT model for 
the two spectra are compatible at 90 % confidence level. 
In order to better constrain the optical depth, r, for the 
simple difussion model (see Discussion) we fixed the seed 
photons temperature kTo and the hot plasma temperature 
kT e to a single value (see Table 1). In Fig. |3we show the 
unfolded EF(E) spectrum of the entire observation, plot¬ 
ted together with the residuals in units of a with respect 
to the best fit COMPTT model. The best fit parameters 
together with the 90% c.l. errors for a single parameter 
are reported in Table 1. 


3. Discussion 

We analyzed for the first time the shape of the pulse pro¬ 
file of SAX J2103.5+4545 at energies above 20 keV. The 
observed X-ray pulse profile shows a broad main peak and 
a secondary peak that seems to become more prominent 
at higher energies. The differences between the spectra of 
the main (pulse phase: 0.5-1) and secondary (pulse phase: 
0-0.5) peak are evident in their ratio, which increases with 
increasing energy up to ~ 30 keV and then decreases, in¬ 
dicating that the main peak is characterized by an excess 
of emission in the energy range 20-50 keV; the difference 
in the mean 5-200 keV luminosity between the two pulse 
phase resolved spectra is a factor of about 1.5. The pulsed 
fraction is observed to significantly and steeply increase 
with energy, from ~ 45% up to 40 keV to ~ 80% in the 
range 40-80 keV. A similar behaviour was also found for 
the Be/X-ray transient RX J0117.6-733 where the pulsed 
percentage increased to at least 79% in the 20-70 keV band 
(Macomb et al, 1999). 

The average spectrum, as well as the pulse phase 
resolved spectra (i.e. main peak spectrum and sec¬ 
ondary peak spectrum) can be well described by a 
Comptonization model. According to this model, the seed 
photons for the Comptonization have a temperature of 
kTg ~ 1.5 keV, injected from the bottom of the accretion 
shock of electron temperature kT e ~ 13 keV and Thomson 
optical depth r ~ 2.4. The electron temperature is too 
high for an accretion disk ( [Frank, King fc Raine 2000| . 

In the hypothesis that the two peaks in the pulse pro¬ 
file come from the two polar caps on the NS surface, the 
variation of the pulse profile with energy cannot be simply 
explained by geometrical effects (the viewing angles of the 
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polar caps, when directed towards the observer, are differ¬ 
ent and the projected emission areas in the line of sight 
are different) because the geometrical effects can change 
the observed flux but not the spectral shape. We there¬ 
fore propose that the differences between the spectra at 
the two peaks is mainly caused by scattering effects. Since 
the main peak do not vary significantly with the energy 
we consider in our model that only the second peak is 
scattered. 

Considering that the radiation in the main peak comes 
from one polar cap and from the continuum of the NS, 
the flux can be written / spot + /ns- The secondary peak 
then is I' svot e~ T + /ns (part of the secondary pole emis¬ 
sion (1 — e~ T ) is scattered away from the line of sight 
and does not contribute to the the observed peak flux), t' 
is given from the Comptonization fit parameter and from 
the Klein-Nishina equation, noting that in our case kT e <C 
me 2 . Above 30 keV, the pulsed fraction is 80% which 
shows that the continuum from the NS has a small contri¬ 
bution. In that energy range, the ratio given in Fig. [^rep¬ 
resents the ratio of the two flux peaks, which indeed is well 
fitted by the simple diffusive model, (1 + a)/(l + a'e~ T ), 
where a = / spo t//NS and a' = /' pot //NS (note that / spot 
and /' pot may be different due to geometrical effects, since 
the viewing angles of the magnetic poles, when directed 
towards the observer, may be different giving a difference 
in the projected emission areas along the line of sight). 
Although the statistics are not good enough to constrain 
all the parameters (the error bars on the single parame¬ 
ters are quite large), interestingly we always obtain that 
a' = 0.1a. We can therefore quantify the geometrical dif¬ 
ference between the two projected emission areas in 10%. 
On the other hand, in the energy range between 5-30 keV 
with a pulsed fraction of ~ 50%, the continuum of the 
NS contributes half of the emission. The emission of the 
secondary peak is likely scattered through the Compton 
shock layer, above the hot spot, and contributes to the 
continuum. 

Finally, from the observed double-peaked folded light 
curve, we infer a limit on the angles that determine the 
geometry of the accretion powered pulsar (if the NS has a 
pure centered dipole field); these are the angle between 
the rotation axis and the direction to the observer, i, 
and the angle between the magnetic and rotation axes, 
(3. Assuming a typical mass value of 20 M 0 for the BOV 
companion star (Reig et al. 2004) and a 1.4 M 0 NS, from 
the mass function the inclination angle of the system is 
i ~ 30°. The parameter space when two polar caps are 
seen is determined by the condition i > 90° — p — Sg, 
where Sg is the gravitational light deflection of the emit¬ 
ted X-ray ( Beloborodov 200 2+ In our case, considering Sg 
= 0, this leads to (3 < 60°, or higher with non zero Sg. 
The gravitational light bending effect is less then 10° (see 
Bulik et al. 2003). 
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